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Abstract 

Two-point correlation functions of spin operators in the minimal models M PiP i 
perturbed by the field $13 are studied in the framework of conformal perturbation 
theory pQ. The first-order corrections for the structure functions are derived ana- 
lytically in terms of gamma functions. Together with the exact vacuum expectation 
values of local operators, this gives the short-distance expansion of the correlation 
functions. 

The long-distance behaviors of these correlation functions in the case M.2,2n+i 
have been worked out using a form-factor bootstrap approach. 

The results of numerical calculations demonstrate that the short- and long- 
distance expansions match at the intermediate distances. Including the descendent 
operators in the OPE drastically improves the convergency region. The combination 
of the two methods thus describes the correlation functions at all length scales with 
good precision. 
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1 Introduction 



A complete set of correlation functions is the main object completely characterizing field 
theories. The well-known nontrivial examples where exact correlation functions have 
been found are the two-dimensional scaling Ising model in the zero magnetic field and 
the conformal field theories PUIS!- Unfortunately, away from the free fermion/conformal 
points, the specific methods that have been applied for those models are difficult to 
generalize. In particular, the problem of finding exact analytic expressions for correlation 
functions of integrable two dimensional massive models is still open. In approaching 
this problem, we follow [I] and study general behaviors of correlation functions in the 
perturbed conformal field theories at all scales by applying a combination of conformal 
perturbation theory and the form-factor bootstrap approach. 

In what follows, we concentrate on the minimal models of CFT M. PtP ' (with coprime 
integers 1 < p < p') perturbed by the field $13. These massive field theories having an 
infinite number of conservation laws are integrable [3]. From the statistical mechanics 
standpoint, the given massive models describe a universality class of the integrable RSOS- 
type models |3] in the corresponding antiferromagnetic regime. 

The conformal field theories M. P)P > underlying the massive models under consideration 
are characterized by the central charge of the Virasoro algebra 

c=l-6 (p '-f . (1) 
pp 

There are (p — 1) x (p' — l)/2 primary fields = 1, • • • ,P — 1 and k = 1, . . . ,p' — 1) 

in the model. The conformal dimensions (A/^, A^) of these fields are determined by the 
Kac formula: 

(p'l -pkf- (pf - pf 
= 4^ ' (2) 

In what follows, we impose the standard normalization for the primary fields adopted in 
CFT, 

($i,*0*0$i,k(°)>CFT = M~ 4Ai ' fc • (3) 

We use the symbol / = $n for the unity operator throughout this work. The operator 
^> = $12 is identified with the spin operator, and $ = $13, with the energy operator. It 
is convenient to use the positive rational parameter 

p — p 
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instead of p and p' because, for example, the operators J, $, and $15, which are 
important for us, then have the conformal dimensions 



i - 2 i - 1 4£ - 2 

The general procedure developed in E] allows computing structure constants in the 
operator algebra and determining all correlation functions. This knowledge can be used 
to study correlation functions of the $13 perturbation of the model M. p>p i as follows. 

The scaling model we study can be formally defined through the action 

A Spp , = A Mp p , +g j ' §d 2 x , (4) 

where Am / denotes the "action" of the critical model and g is a coupling constant having 

the dimension g ~ m^~ 2A *, where ni\ is a mass gap (the mass of lightest particle). We 
note that the operator $ has the conformal dimension A$ < 1 and the perturbation is 
relevant. More explicitly, the dimensional parameters g and mi are connected via the 
exact "mass-coupling-constant relation" found in jH]: 

*9 = ~ , ^t^ 2 , U—M—j)*rn& . (5) 
y (£- l)(2f- 1) \ K Z + 1 J + 

Here and hereafter, the massive parameter m is related to the mass of the lightest particle 
rrii by 

'*7(i - fh(¥) 
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m = mi 
and we set 



70) 



8 sin 7r£ 

r(x) 



Formal definition (@J) is understood in the sense of the perturbation series. For exam- 
ple, for the correlation function of the spin fields, we have 

(*(*)*(0)> = / (^(x)^m(yi) ■ ■ ■ Hyn))cFTd 2 yi ■ ■ ■ d 2 y n . 

n=0 ' J 

Here the symbols ^ in the l.h.s. denote the scaling fields in the perturbed model; in the 
UV limit, these fields become the conformal spin fields $12. Throughout this work, we 
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generally assume that the renormalized fields in (jlj) have definite scaling dimensions and 
are denoted by the same letter as in the conformal g — > limit. We also require that the 
normalization of scaling fields at the limit be fixed as in 

The UV and IR regularization schemes for the integrals over the plane were discussed 
extensively in pQ. We only briefly recall the final prescription here. Since the perturbing 
field is relevant, the UV renormalization can be achieved by adding a finite number of 
counterterms. To handle the IR divergences and to work out the infrared-safe perturba- 
tion theory, the so-called conformal perturbation theory is developed pQ. 1 The conformal 
perturbation theory is based on the following hypothesis. Let {A a , a = 0, 1, . . .} be a 
complete set of local fields in the scaling theory. It is assumed that if the renormalized 
fields that are eigenvectors of the dilatation operator, i.e., the fields with a definite scaling 
dimensions, are chosen as basis elements, then the structure functions in the correspond- 
ing operator-product expansion 

A m (x)A a (0) = J2ctA a (r)A b (0), (6) 

b 

are analytic functions of the coupling constant g. This is a rather natural conjecture since 
the functions C A h mAa {r) as local quantities are assumed to not develop any nonanalyticity. 
We propose that under this choice, the renormalized fields Ab(r) turn out to be pertur- 
bations of the corresponding basis fields (primaries and descendants) of the conformal 
model. Then the functions C A A (r) are given by regular expansions 

Ct Aa (r) = r 2 ^-^-^) (cZ Aa + gr>- 2A QM + (gr^fQ® + ...). ( 7 ) 

The zeroth-order terms here are determined by the structure constants C^ b A from the 
critical model, 2 and the problem of perturbation theory is to define (^-independent cor- 
rections Q®. 

To define the correlation functions following (JUJ), we must also find vacuum expectation 
values (VEVs) of the local fields A a , which might be nonzero in the off-critical case. These 
important quantities are nonlocal and, generally, nonanalytic in the coupling constant 
g. From counting dimensions, we have (A^) ~ m 2Ab . The set of VEVs contain all 
nonperturbative information on the theory. The VEVs for perturbed primaries were 
proposed in ^H] (also see ^H|). The Lukyanov and Zamolodchikov formula for the VEVs 
in the first ground state [5] can be written as 

(^=(-1)^(1 - - 1)) • (8) 

^ee |13| for the conformal perturbation theory in all orders. 

2 We note that because of fusion rules, many of the C A b A in the conformal case are zeroes. But some 
of the corresponding C A b A (r) may become nonzero because of the correction terms. 
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The function Q{rf) is defined in terms of the following exponent of an integral that is 
understood in the sense of analytic continuation from the region where it converges: 



Q(rj) = exp 



dt /cosh(2t) sinh((?7 
TV 



l)t) sinh((?7 + l)t) 



2 cosh(t) sinh(£t) sinh((£ + l)t) 



77 2 - 1 



-2(C+l)t 



(9) 



We note that if £ > 1, for instance, as in the models M. PtP+ i of principal series perturbed 
by the energy operator, then there are no bound states. In that case, we use the kink 
mass in the expressions like (JKJ). 

The short distance expansion for the two-point correlation function of spin fields is 
determined by the VEVs of fields $ifc (and their descendants) with odd k = 2s + 1 (see 
Eq. (|T3j) below). After some computations we found that for these fields the expression 
simplifies and can be written explicitly in terms of gamma functions as following 



Q(i-2 S o=nQi(o 



(10) 



where 

Qi(0 - 

x 7 



[l- {21 - 1)0(1 -2ZQ w 

V ^> 

2(£ + l) 

(2l-l)£\ fl-(2Z-l)£ 

7 



7 



l-2(/-l)e 



7 i- 



(2/ + i)e 



(ii) 



The one-point VEVs for the first nontrivial descendent operators are also known in 
the analytic form |2()j . 



(L_ 2 L_ 2 $ lfc ) 
where the function W(r/) is 

1 



-(l+0 4 W(l-£(£:-l))m 4 (<l> lfc >, 



(12) 



W(r/) := 



£ 2 (£ + l) ; 



7( o )7(— J7( o h( R— > 



Although the problem of determining exact VEVs of all descendants is still open, the 
knowledge achieved up to now can already be used to analyze correlation functions. For 
example, VEVs (|8 |) -([12 p determine the leading terms in the short- distance expansion of 
the two-point correlation function of spin operators, 



(*(*)¥(0)> = Ck(r)<I) 



^(r)($(0)) + C*-(r)($ 15 ) 



'\J>\J> I 



r)(L_ 2 L_ 2 /> + C^ 2£ - 2 *(r)(L_ 2 L_ 2 $(0)> + 



(13) 
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which is the main object of study in the first part of this work. The conformal pertur- 
bation theory based on an exact knowledge of VEVs of local operators is thus applicable 
for an effective study of the short-distance (mr <C 1) behaviors of correlation functions. 

On the other hand, theories (@J) are massive integrable models with factorized scat- 
tering. The scattering matrices in those theories coincide with those in the restricted 
sine-Gordon model [TQl HH E2] • We can apply the spectral expansion |HJ E] (see [T7] and 
the references therein for details) 

<*(*)»«>)> = £ £ k I fr • • • fr<°i»<°>iA' • • ••*>.-» x 

n=0 {a„} ' J 

x 01 ..^ (/?!,..., AO|¥(0)|0) e-^^ cosh ^ (14) 

for the correlation functions and study the same correlation functions at long distances 
mr ^> 1. Here, are the particle types, and the quantities 

(0|*IA WU. ai ... an (f3 u . . . , (3 n )\m 

are matrix elements of the local operators in the basis of asymptotic states, the form 
factors. The latter are assumed to be meromorphic functions of the rapidity variables f3 
satisfying a set of locality axioms (7j. For the case of primary fields, the exact solutions 
for the form-factor equations have already been found [TJ [TUl 1211 EH 121] • 

The difficulty not easily handled in the form-factor approach is that the axioms for 
scalar primary fields do not fix an overall normalization of the form factors or, in other 
words, the normalization of local operators. To impose conformal normalization (J3J), 
we determine the overall rapidity-independent factor as in QjJ |2H] by the VEVs of local 
operators (JBJ) (also see |30j). This prescription seems quite natural from the standpoint 
of the quantum group procedure |%1 ITUl IP3]. 

In the integrable cases, there are thus two different expansions (|13)) and (|14j) of the 
same correlation functions of the local operators. A natural idea would be to compare 
these formulae for short and long distances. As soon as the expansions are matched 
somewhere at the intermediate distances, we would obtain a description of correlation 
functions at all scales. The practical question, of course, is how easily that can be 
achieved. 

We mention that this program of studying correlation functions, initiated in (also 
sec [20 ) for the case of the Lee- Yang model, has already been applied to several more 
complicated cases [321 EH EH] • 

We have found that for the models M PjP ' perturbed by the field $ 13 in the resonance- 
free cases, an elegant analytic form for the first-order perturbation theory for the structure 
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Figure 1: Correlation function of the two spin operators (in units mr) in the minimal 
model M.2,7 perturbed by the field $13. 

functions can be easily provided. With answers (jffi-(|12|) taken into account, this leads 
to expression controlling the UV behaviors of the spin correlation functions. 

For the case M.2,211+1, this together with the form-factor expansion up to three particles 
already gives a broad region where the two expansions match. 3 For example, one of the 
plots illustrating this phenomenon for the model M.2,7 perturbed by the energy operator 
is shown in Fig. ^ The large region of convergence between 1 and 3 makes it clear that 
the given method provides a good approximation for correlation functions at all scales. 

This matching between different expansions at the very first steps demonstrates that 
the combination of the conformal-perturbation and form-factor approaches is a very useful 
tool for checking the consistency of different proposals used in both constructions. For 
instance, in our case, we obtain an additional confirmation for statements that the form 
factors for the primary fields are given by the minimal solutions with the prescribed 
analytic properties [2*7| I29j. that the normalization of the form factors is fixed correctly, 
and so on. Moreover, that the contribution from the descendent operators drastically 

3 We choose the model M.2,2n+i just for simplicity. The form factors of primaries in this case do 
not contain contour integrals. Another perturbed CFT with that property is a parafermionic model 
perturbed by the first energy field 
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improves the convergence of the two expansions (see Fig. ^) can be treated as a good 
check for the correctness of (|12|). In addition, this method can also be treated as an 
alternative way for studying universal scaling behaviors of the correlation functions in 
the corresponding lattice RSOS models [31 153]. 

Although many technical components that we used are already known, we collect all 
necessary facts for the reader's convenience and to make the text self-consistent. This 
work is organized as follows. In Sec. 2, we consider the procedure for finding first-order 
corrections to the structure functions. We demonstrate that the double integrals over a 
complex plane in the physically interesting situations can be taken analytically leading 
to an elegant answer in terms of gamma functions. We thus obtain the UV expansion 
for the correlation function of spin operators. We provide some details of the numerical 
computation for M.2,7 as an example. In Sec. 3, we recall the scattering data and form- 
factor axioms applicable to our case with many particles. We restrict our attention to 
the three-particle case, which already gives good IR expansion formulae. In Sec. 4, we 
give some results of the numerical computations and demonstrate the consistency of both 
expansions. In Sec. 5, we draw final conclusions. The technical details are collected in 
the appendix. 

2 Short-distance expansion scaling M PI / model 

Our main aim is to study two-point correlation functions of the spin operators in the 
massive $13 perturbation of M. PtP ' model (jU). In this section, we analyze the short- 
distance expansion based on the operator algebra 

(*(aO*(0)> = CUr)(I) + Cf*(r)<<&(0)> + C*Jf (r)($ 1B > 

+ C^ L - 2 \r){L. 2 L^I) + C^ 2£ - 2 *(r)(L_ 2 Z_ 2 $(0)> + • • • . (15) 

In the r.h.s., we keep only the leading terms in the small r that have nonzero VEVs. We 
thus omit the derivative operators as well as operators with nonzero spins. Taking the 
known VEVs of local operators (JSj)- (fT2j) . we study correlation functions (fT3|) following 
the program initiated in pQ. 

2.1 Perturbation theory for structure functions 

According to PQ, the structure functions in (|15|) have the expansion with the zeroth order 
given by the structure constants from the underlying M. v ^ (J7J), 

C&(r) = C*, r 2A -~ 4A * + Cf„ (1) (r) + ■ ■ ■ . (16) 
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The first-order corrections (r) can be represented in the regularized form with R 

being the parameter of the infrared cutoff, 

C^ (1 \r) := lim \-g [ (A K (ooMy)^(x)^(0)} CFT d 2 y 

«->°°L J\y\<R 

+7ro C^Cqa ^2(A K -A A ~A $ +l) r 2A A -4A^l 

~* A K — A A — A$ + 1 J ' 

As in the Lee- Yang case PQ, it is convenient for practical purposes to apply an analytic 
regularization procedure for the integral. In our case of two operators \I/ = $12, we 
analyze the regularized integrals 

C$i\r) = -g J (<f>(yMxW0)) CFT d 2 y, 

C^\r) = -9 f Woo)*(yMxmo))orr*y > ( 18 ) 
C|^ 5(1) (r) = -g y"($ 15 (oo)$(y)*(^)*(0)) CFT A . 

We recall the exact form of the conformal correlators appearing here. The correlation 
functions in the integral in the first line above are defined via the structure constants of 
the operator algebra 

($(y)*(x)¥(0)) CPT = C^\y\- 2A *\x\- 2( - 2A *- A ^\y - x|" 2A * . (19) 

The correlation function including the field $15 has a very similar form. Indeed, the 
correction term in the third line can be rewritten as 

C% 5(1 \r) = _^ r 2(l-A s -2A, + A 15 ) 



/ 

2- 



x / |z| 2(A * +2A *- 2 - Al5) ($ 15 (oo)$(l)^(z)^(0)) CFT ^ 



where 

($ 15 (oo)$(l)^(z)M/(0)) CFT = C* 9 C*$ \z\& |1 - z\& . (20) 

In general, the four-point correlation functions of the primary fields in the conformal 
models can be found, for example, using the free-field construction of Dotsenko and 
Fateev [H] or the differential equations method. For example, the integral in the second 
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line of Eq. (|18|) is expressed through the correlation function 

(tf(0)^(z)$(l)$(oo))cFT 



i ; 'U + l' £+1 e+1 

2 2£ 2 



z A(i+T) (1_ 2;) 2(«+i) i? 7—^, T 

With Eq. (JZEJ), it has the form 

Cji 1} (r) = -^r 2 ~ 4A * y" | 2 | 2A *- 2 (vI/(0)vI/(^)<l>(l)$(oo))cFTrf^ 



«±4 717+t) ,. ( f l2 /" 12 I l-^M 1-^- 

7 2 (|+T)-°vy y 



X M 5 |1 — 2/p + F — 2/ 



2S 



2c 1 



2S 



(21) 



(22) 



Here, we fixed the way of an analytic regularization in the second line that allows com- 
puting the integral explicitly before taking the limit. To do this, we use formula (j2HJ), 
which is derived in the next section. 



2.2 The integrals 

From explicit formulae (|19j) - (j21j) . we find that the properly regularized corrections (|18j) 
are reduced to integrals of the forms 



JiM = j d 2 y \y\ 2p \y-M 2 \ (23) 

J 2 {a,b,d,e,c) = [ d 2 x [ d 2 y \x\ 2a \l - x\ 2b \y\ 2d \l - y\ 2e \x - y\ 2c . (24) 



The first integral can be easily computed by reducing it to a product of beta-integrals, 
giving 

) = H (^M£±i). (25) 

7(p + q + 2) 

Integral (}2"4*|) is more complicated. It can be rewritten ^SJ E3 121] i n terms of a prod- 
uct of two contour integrals admitting a representation via the higher hypergeometric 
functions 3F2 at unity [25J (also see [34J). For completeness, we collect the details of 
computations in Appendix B. 
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In our cases (J17|) . the integral J2{a, b, d, e, c) appears at special values of the parameters 
and can be expressed in terms of gamma functions due to the following formula 

f p-a-1 5 /3-a-l 

J 2 , o — p — 1, a — 1, 



X7 (^|±l) 7 (£z£z^±l) . (26) 

Using (|2"5|) - (j26j) . we can analytically derive first corrections for the structure functions 
for arbitrary £ > 0. 

2.3 First-order corrections to the structure functions 

Applying the formulae for integrals (|25 p -([26 p and the exact values of the CFT structure 
constants from Appendix A, we now give exact expressions for the structure functions up 
to the first order in g. 

The functions C^(r) and C*^ 5 (r) are the simplest. (We recall that the conformal 
field theory structure constants were = 1 and = 0.) Using ®, we find that 
they have the forms 



V7(#i)7(g)7 2 (^) 



C$? (r) = -ng r W C* 9 C*$ . (27) 

The first-order correction to the function C*^(r) 

C^'{r) = -grt+i \ hmJ(e), 
7 2 (^t) ^° 

t/ \ w £ 2£ £ + 3 £ £ 

./(e) = J 2 (-— — + 2e,——-e, 



is computed by applying Eq. (|26|). Taking the limit, we obtain the expression 

w v 2(e + i) 2 ^ 
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Finally, the analytic expression for the structure function C*^(r) in the first-order per- 
turbation theory is 



" () ~ U(£h(f3)J + 9 ( -WTW) ' (28) 

It is convenient to rewrite structure functions (|27j) - (|28jl in terms of the lightest particle 
mass. Applying the mass-coupling-constant relation, we have 

Ct.(r)=rBJl- (mr )A « + 1 > 2 M^WM^ YI 
1 1 ' (f-l)(2«-l)U(^Wt^MjTl) 



.t(^t)t(H 




est,) = - rW(mr) A _ i|_|L _ (^ ^p^ ) ■ (3D 

We also need the structure functions for the descendent fields in the main order. These 
are expressed using the corresponding CFT structure constants of the primary fields just 
as 

4(£ + 3) 5 



C^ L - 2l (r) = 77 * r -W-C74» , (32) 



2.4 UV expansion for the correlation function of spin operators 

Taking these formulae for structure functions, VEVs (JSJ) as well as the expressions 

(L_ 2 I_ 2 /) = -m 4 (^-) 2 , 
12 



for VEVs (|12|) into account, we obtain the explicit first-order UV expansion in g for the 
correlation function 

(tf (r)tf (0)) = rfrr [A(r) - B(r)Q(l - 2f) + D(r)Q(l - 40] . (35) 

We recall that the function Q(rj) © for the fields $i,2 S +i is given by the Eqn. (fTUj) - (JTTJ. 
For the fields $13 and $15 the respective expressions in ()35|) turn out to be 



(1-0(1-20(1-30(1-40 



4(e + l) 2 



X7 



7 



1 - 2£ 



7 



l-3£ 



7 



e+i 



7 



1-4^ 



The r-dependent functions A(r), B(r), and Z}(r) in ()35|) have the forms 





A(r) 



7 2 d: 



mr 



4(^ + 3) 2 7 2(i±i) 

(e + i) 2 M^b^Mgr) ' 

(1-0(1-20 I 7(^)7(gf )7 2 (^t) 



mr w 



B(r) 



1 + 



(e + 1) 2 7(1)7^) 
(e-l) 2 (3C + l) 2 7(l)7(¥) 



'7 8 (^)7 9 (^)7(^t; 



D(r) 



7(^)7(rfi) 

J(ir)7@), 
0(i-Q 3 

4(e+ i) 2 (i-20 

0(i-0 

4(l-20(3e + l) 2 U(£) 7 (^) 7 (*=S) 



[mr] 
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2^2-2^ 
€+1 ' 



7 



7/ -J-Wf kJ. 



2(1-1) 

lmrH+ 1 > (mr) £+ J 



8| 

fmr)^ 1 . 



(37) 



(38) 



(39) 
(40) 



2.5 Model A^j as an example 

The first model in the series A^.s is identified with the Lee- Yang model ^1] with the 
central charge of Virasoro algebra c = — y. There are only two primary fields in the 
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theory: the unit field / := $xi — ^14 with the conformal dimension zero and the spin 
field $12 = $13 with the dimension — ~. This example was analyzed in detail in PI 1201 • 
The next example is the model M.2,7 [15 j with the central charge 

c=-f. (41) 

In this case, there are three primaries: the unit field / := $n = $16, the spin field 
^> := $12 = $15, and the energy field $ := $13 = $14 with the respective conformal 
dimensions 

A/ = 0, A v = ~, A* = -j. (42) 
The first-order corrections in g to the structure functions are given by expressions (|2T)|) - 

J2D with £ = §. 

Providing an analytic continuation of the Lukyanov-Zamolodchikov solution, we have 
the following data for the VEVs of operators giving the leading contribution: 

mf ($13) = 2.269550689, 
m f ($ 15 ) = -2.325136 i 

for primary fields and 

mfc^ 2L - 2l (r)(L_ 2 l_ 2 I) = 



-0.0005899419474(mir 
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30 

= 0.01123584810(mir)- 

for the first nontrivial descendants. Taking into account that the mass coupling constant 
for f = 2/5 is 

g = -0.04053795542378225mf /7 (43) 
and collecting all data together, we obtain the short-distance expansion up to r 6 : 

m\ (*(ar)*(0)) 5.83(mir)? + (mir)? - 0.265(mir) 22/7 + 0.233(m!r) 24/7 

-0.045(m!r) 4 + 0.011 

-0.0006(mir)^ + 0(r 6 ) . (44) 

The plots for this function are shown in Figs, It can be clearly seen that the 

contributions from the descendent fields (L_ 2 L_ 2 Q) and (L_ 2 L__ 2 I) are essential. 
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3 Long-distance expansion 



Massive models (@J) allow an alternative description as scattering theories with factorized 
S matrices (TU1 HI] . In this section, we use the form-factor approach to study the long- 
distance behaviors of the spin correlation functions. The structure of ground states of 
the model as well as the spectrum of particles depend on the numbers (p,p') CHI H2j • 
There are two kinks interpolating between p — 1 different vacuums and [(p' — p)/p] bound 
states. 

For simplicity, we restrict our attention to the case of nonunitary models M.2,2n+i 
having only n primary fields $1^ = 3>i,2n-fc+i> (k = l,...,n). The corresponding per- 
turbed models are the simplest integrable two-dimensional models in the sense that the 
scattering matrices are diagonal. The multiparticle form factors of the fields \&, $, . . . 
in these theories do not contain complicated contour integrals originating from solitonic 
form factors jjj- Hereafter in this section, we use the symbol £ only for 

e = -H_. 

S In - 1 

We note that for studying a general case, one must also take kink contributions into 
account. 



3.1 Perturbed A^2,2n+i models as scattering theories. 

We briefly recall the basic facts about the scattering theories under consideration (see, 
e.g, |27 a for a detailed review). In A^2,2n+i models perturbed by the field $, there are 
j — 1, . . . , n — 1 scalar self-conjugate particles with the masses 

sin(^) , 

m i :=m i ■ 7^ • ( 45 ) 
sm(f ) 

Here, m,\ is a mass of the lightest particle "1". The standard momentum parameterization 
is given in terms of rapidities j3 as 

p° a = m a cosh (3 , 

p\ = m a sinh (3 . (46) 

The S matrix for these models originates from the breathers S matrix of the corre- 
sponding sine-Gordon model [TlH 12*6]: 

min(ab) — 1 

S Qb (J3) = f ] ^ L (j3)f ! ^(J3) J] (/i°W/?)) 2 - (47) 

s=l 
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Here, the basic building block f a (f3) is defined through the meromorphic functions 

tanh^ + ^) 
tanh ~(p — iira£) 

In the physical strip < lm/3 < ir, the scattering matrices satisfy the unitarity and 
crossing symmetry conditions given respectively as 

S ab ((3)S ba (-P) = 1 , S ab ((3) = S ab (i7T - (3) , 

as well as the additional bound-state condition described as follows. If S a b(l3) has a simple 
pole j3 a t, = iu c ab in the direct channel, then the particle "c" of mass 

m c = m l + m t + 2m a m b cos(u c ab ) 

is a bound state of "a" and "b". The bound state \B C {(3)) with (3 = [3 a + (3 b is defined as 
the projection on the pole of the two-particle state \B a (p a )B b (P b )) at the relative rapidity 
Pab = iUab- Then the S matrices for the composite particle can be found via the bootstrap 
equations 

S cd (P) = S ad (P + i(n - u b J)S bd ((3 - i{n - <)) . 

In particular, the total scattering matrix S ab (jTrj) can be reconstructed starting from the 
S matrix Sii(/3) = for the lightest particle assuming that the bootstrap tree for our 
massive models is closed under the fusions [2E] 

cii x dj — > a>i + j or di x aj — > a-2n-i-i-j , 

The rule in the first line above is to choose the variant where a final particle i + j or 
In — 1 — i — j is in the range (1, . . . , n — 1). 

3.2 Form-factor approach 

As we already mentioned, the correlation functions admit a form-factor expansion that 
is very useful for studying the infrared behaviors. For example, for spin fields, the corre- 



16 



sponding spectral decomposition is 

(*(*)*(°)> / ^<°I*(*)IA, • • • , Pn) ai ... an X 

n=0 ' ^ 

x ai _ an (p 1 ,...,p n \^(0)\0) = 

EJ_ /" ^ ^-rE.-w.; cosh /3j 

„.| Z-/ / 9-tt " ' 



nl^J 2tt 2tt 

n=0 {aj} 



xF a „... ai (/?„,..., (/?!,..., n ) , (49) 

where we introduce the matrix elements of the local operators in the basis of the asymp- 
totic states 

F ai ... an (f3 u ...,&) = (0|¥(0)|fr, . . • , /3 n ) ai ... 0n . (50) 

Form factors (|5U|) are defined in this approach as a set of functions satisfying Smirnov's 
axioms [7] (also see jH])- We recall them for completeness. The first two axioms are the 
Watson equations 

F ai ... ajdj+x... \P%i • • • j Pi: Pi+1: ' ' ' i Pn) 

= Saja,j-\-\\Pi (3j+l)F ai ... a . +iaj ... an ( y (3i, . . . , Pj+i, Pj, ■ ■ ■ , Pn) i 
F ai a 2 -a n (Pl + 2ni } p 2} ■ ■ ■ ,P n ) = 

= F a2 ... anai (p 2 ,--- ,p n ,p 1 ) . (51) 

The requirement for relativistic invariance implies that for the local operator with spin 
s, one has 

F ai ... an {Pi + A, • • • , p n + A) = e sK F ai ... an {Pu ••-,&)• (52) 

The ultraviolet bound condition requires that there exist a finite constant t(J, n) < oo 
such that 

F ai -a n (Pi + A, • ■ ■ ,Pj + A,P j+1 ,...,p n ) = 0(expe' ( ^ )|A| ) , A — ► oo . (53) 

In addition, there are two pole conditions. The first is a kinematic pole condition. In our 
case for self-conjugate particles, it looks like 

-i lim (^' - P)F aaai _ an (p' +i7r,P,Pi,...,P n ) = 

n 

= (1 - ]JS aaj (P - P,))F ai ... an (p u ...,Pi). (54) 
i=2 
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The second pole condition is the equation for bound states. Let <2j x aj 



Ofe. Then 




(55) 



where the three particle on-shell vertex is defined as 4 



iRes S ab ((3) 



(K b ) 2 ■ 



ab 



Smirnov [7] proved that the set of functions (J5U|) satisfying conditions (jB*T|) - (j53|) defines 
a matrix element of a local operator in the scattering theory with the given S matrix. 

3.3 First form factors 

The solutions of locality axioms (joT|) -(f53 j) for the primaries in perturbed M.2,2n+i mod- 
els have already been discussed in the literature. The form factors of the fields \1/ and 
$ were found in ^U] via the quantum group reduction procedure from the sine-Gordon 
model. The case of the scaling Lee- Yang model M.2,5 was directly elaborated in [T]. The 
general expressions for the primaries $1^ appeared in j2Z] based on the form factors for 
exponential operators in the sinh-Gordon model f2Bl. But for our aim of studying the 
correlation functions that are in agreement with (j!5p . we must complete these prescrip- 
tions by choosing a proper overall factor determining normalization Q of local operators. 
This was essentially provided in (2H|. 5 Indeed, after the quantum group reduction |51IT3]. 
the expressions for the first breather form factors of the exponential operators in the 
sine-Gordon model can be essentially treated as the following particle "1" form factors 

4 We note that the perturbed minimal models A^^n+i violate the one-particle unitarity since some 
of its couplings arc purely imaginary. 

D Another advantage of Lukyanov's bosonization rules is that they can be easily extended to include 
the kink sector for a perturbed minimal CFT with general 



18 



of the primaries in the reduced model: 



.6 



(0|$ lfc |0) = ($ lfc ), (56) 



sin ((k - l)</2 



(0|$ lfc |/3) = i C v . ^ ($ lfe ) , (57) 

sm(7r£) 



sin 2 ^(A; - 1)tt£/2 



(0|$ lfc |/3 2 ,A) =i 2 C 2 \ 2 ^ R(Pi-fo) (*ik), (5 

sin (7I"£J 

sin f(jfe - 1)^/2) 

(0|$ lfe |/3 3 ,/5 2 ,/3i) = z 3 C 3 V . ^ TT^.-ft) 

sin(7rt) - LJ - 



sin 2 f(Jfc - l)</2 



sm 2 K) ns< .2cosh(^ 



<*ife> , (59) 



and so on. Here, the function R(/3) determining the rapidity-dependent part of the two- 
particle form factor is given explicitly as 

r,/^ f „ /"°° dt sinht sinh£tsinh(£ + l)t , . z' , | . „, 

R{/3)=expU .\ 2o / L cosh2 (60) 

We note that it has a property that is very useful for checking the locality axioms at low 
particle levels (while one of the easiest ways to prove the validity of the general formulae 
is to use the free-field realization): 

WW***) 



sinh(/3) =F i sin(7r£) 



n n (, + M (a + ,) R (, _ Ei (a + 1} ) = . (.i> 

The /3-independent constant C fixing the normalization of the first particles is 

C 2 = 8cos 2 (^)sin(^)exp(- f* ~) ■ (62) 
2 2 V y 7r smt/ 



6 We note that the same result can be obtained using the exact realization of the form factors in the 
RSOS model and taking the scaling limit 
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3.4 Correlation functions at long distances 



One can see that the leading terms for the correlation function already come from the 
zero- and one-particle form factors. For example, even the first two-term expression 

($ifc(x)$i fe (0)) = 

2 r 2 sin 2 (|(fc-l)Q / r<dt t \ X 

= ($i fc il — — exp - / — )K (m l r)\ + ••• , 63 

I 7r sm(~£J V Jq 7T suit/ J 

where K (x) is a Macdonald function, turns out to be a good approximation for the 
correlation functions at long distances. 

Of course, for more accurate calculations, one must take higher form factors with 
more particles into account, including those with particles a,j with j > 1. The latter can 
be easily derived directly through the bound-state axiom (|HK|) . 7 It is clear from (|4*9*|) that 
at long distances the contribution of an n-particle form factor with particles 
and masses m ai , . . . , m an to the correlation function is essentially approximated by the 
quantity exp(— r Y^j ^cy)- In practice, the form factors with a higher number of particles 
ai and particles of high mass become important only at small distances. It is useful 
to order the form factors according to the value of this quantity. Ror example, for the 
perturbed M.2,7 model that is considered further, the leading terms thus come from F\, F2, 
Fn, F21, Fm, etc., in agreement with mass relations (J43)) mi < m 2 < 2mi < mi + m^ < 
3mi, etc. 



4 Numerical results 

The form-factor and conformal perturbation expansions work well at long and short dis- 
tances respectively. In this section, we collect the numerical data demonstrating matching 
between these two at mr ~ 1 for the case of scaling Ai2j models as an example. Substi- 
tuting £ = I for the M.2,7 case, we obtain the expression 

($i 2 |Bi09))(fli(/?)|$ia) = -($i2> 2 • 0.7574894945 . 
Taking into account the numerical value for the VEV 

ml ($ 12 ) = -2.325598436 i , (64) 

7 The additional requirements for the consistency of formulae H56(l - H59|l . etc., appearing from the 
bound-state condition are automatically satisfied for the operators <&i.fc as was argued in [57]. For form 
factors with a small number of particles, this can be verified straightforwardly. 
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we find that the leading term in the long-distance expansion is given by the formula 

m\ ($12 (a?) $12(0)) 5.408408086 (l - 0.2411163946 • K^rm^ . (65) 

Comparison of the short- and long-distance expansions for the two-point correlation 
function (\l/(x)\l/(0)) is shown in Figs. EtD 




0.00 0.05 0.10 0.15 0.20 



mi- 



Figure 2: Form factor contributions at small distances 

The numerically calculated data is presented in Table The difference between 
the short- and long-distance expansions indicates the self-consistency of the zeroth-order 
perturbation theory and the form-factor expansion up to two particles. 



5 Conclusion 

We have demonstrated that there is a region in the mass scale where the long- and 
short-distance expansions for the correlation functions of the spin operators match each 
other. We can therefore conclude that the combination of the two approaches allows 
reconstructing the correlation functions at all scales. 
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Short 


-distance expansion 


Long- 


distance expansion 


mr 


up to r 8 / 7 


up 


to r 28/7 


up to r 36//7 


up to Fn 


up to Fi 2 


up to Fin 


0.02 


-1.89536 


-1 


89536 


-1.89536 


-1.78911 


-2.14822 


-1.97677 


0.04 


-2.29916 


-2 


29917 


-2.29917 


-2.20386 


-2.41955 


-2.34353 


0.06 


-2.56977 


-2 


56979 


-2.56979 


-2.48822 


-2.63930 


-2.59599 


0.08 


-2.77773 


-2 


77778 


-2.77778 


-2.70767 


-2.82117 


-2.79351 


0.10 


-2.94806 


-2 


94817 


-2.94817 


-2.88736 


-2.97622 


-2.95732 


0.12 


-3.09288 


-3 


09306 


-3.09306 


-3.03986 


-3.11138 


-3.09787 


0.14 


-3.21905 


-3 


21935 


-3.21934 


-3.17245 


-3.23117 


-3.22120 


0.16 


-3.33092 


-3 


33135 


-3.33135 


-3.28975 


-3.33870 


-3.33115 


0.18 


-3.43139 


-3 


43200 


-3.43199 


-3.39489 


-3.43619 


-3.43036 


0.20 


-3.52254 


-3 


52336 


-3.52335 


-3.49009 


-3.52528 


-3.52071 


0.22 


-3.60588 


-3 


60696 


-3.60694 


-3.57702 


-3.60725 


-3.60361 


0.24 


-3.68258 


-3 


68397 


-3.68395 


-3.65692 


-3.68306 


-3.68014 


0.26 


-3.75354 


-3 


75530 


-3.75526 


-3.73077 


-3.75352 


-3.75114 


0.28 


-3.81950 


-3 


82167 


-3.82162 


-3.79937 


-3.81926 


-3.81732 


0.30 


-3.88103 


-3 


88367 


-3.88361 


-3.86334 


-3.88081 


-3.87921 


0.35 


-4.01850 


-4 


02259 


-4.02247 


-4.00626 


-4.01911 


-4.01810 


0.40 


-4.13682 


-4 


14280 


-4.14259 


-4.12947 


-4.13911 


-4.13845 


0.50 


-4.33032 


-4 


34153 


-4.34097 


-4.33216 


-4.33782 


-4.33752 


0.70 


-4.60061 


-4 


62933 


-4.62699 


-4.62274 


-4.62492 


-4.62484 


0.90 


-4.77130 


-4 


82907 


-4.82226 


-4.82018 


-4.82110 


-4.82108 


1.00 


-4.83076 


-4 


90821 


-4.89757 


-4.89617 


-4.89679 


-4.89677 


1.20 


-4.91088 


-5 


03975 


-5.01671 


-5.01633 


-5.01661 


-5.01661 


1.40 


-4.95019 


-5 


14911 


-5.10491 


-5.10540 


-5.10554 


-5.10554 


1.60 


-4.95765 


-5 


24871 


-5.17110 


-5.17252 


-5.17258 


-5.17258 


1.80 


-4.93933 


-5 


34862 


-5.22123 


-5.22367 


-5.22371 


-5.22371 


2.00 


-4.89960 


-5 


45779 


-5.25949 


-5.26302 


-5.26304 


-5.26304 


2.20 


-4.84170 


-5 


58470 


-5.28902 


-5.29349 


-5.29349 


-5.29349 


2.40 


-4.76811 


-5 


73780 


-5.31231 


-5.31721 


-5.31721 


-5.31721 


2.60 


-4.68080 


-5 


92578 


-5.33150 


-5.33577 


-5.33577 


-5.33577 


2.80 


-4.58133 


-6 


15775 


-5.34855 


-5.35034 


-5.35034 


-5.35034 


3.00 


-4.47099 


-6 


44339 


-5.36541 


-5.36181 


-5.36181 


-5.36181 


3.20 


-4.35086 


-6 


79303 


-5.38416 


-5.37087 


-5.37087 


-5.37087 


3.40 


-4.22182 


-7 


21775 


-5.40703 


-5.37805 


-5.37805 


-5.37805 


3.60 


-4.08465 


-7 


72941 


-5.43656 


-5.38374 


-5.38374 


-5.38374 



Table 1: 
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-4.0 
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Long-distance expansion 

Short: up to r 8/7 

cu t t 22/7 
Short: up to r 

cu w t 24/7 
Short: up to r 

cu , . 28/7 
Short: up to r 

Short: up to r 30/7 

Short: up to r 36/7 
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Figure 3: Behaviors of different terms of the short distance expansion 



This also serves as a good consistency check for the correctness of both methods as 
well as the VEVs of local operators. 
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When the paper was completed we learned from V. Fateev that the quantities J-z{a, b, 
d, e, c) in Eqn. (|2*4*j) at the special values of parameters as in the l.h.s. of Eq. (J2*fjjl can 
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be transformed to the class of integrals studied in [3Bj, where analytical expressions in 
terms of gamma functions had already been proposed. We are grateful to V. Fateev for 
informing us about this. 



Appendix 

A Structure constants and conformal blocks 

For completeness, we here give a derivation of structure constants for the subalgebra of 
conformal fields $ife and relevant conformal blocks that are necessary for building the 
corresponding correlation functions in conformal field theory. To reproduce these known 
results 0, we use the method of functional equations [36j. 

We consider the four-point conformal block with the chiral $12 operator 

G(z) = ($ 12 (z)$ Ai (0)$ A2 (1)$a 3 (oo)> , (66) 

where $a(z) denotes chiral parts of the primary field (z,z) (k = 1, ...,p' — 1) with 
the conformal dimensions Ai )7l .. This holomorphic function satisfies the null vector equa- 
tion 

with /i = 2 (2Ai2+i) • ^ e soai ti° ns °f the equation and their monodromy properties can 
be easily found by noting that it turns out to be the Riemann equation 

d 2 u fl — a — a' 1 — 7 — y\du 



d 2 z ^ V z ^ z — 1 / dz 



+(— + ,(,-i) r = ' m 



a + a' + p + p' + -f + y = l 

with the parameters 





Tlx — 1 


a = 


2 P ' 


P = 


2-n 3 


2 P ' 




n 2 - 1 


7 = 


2 P ' 



/ 1 ni + 1 
« = 1 g — P ' 

713 + 2 
7=1 = — P- 
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Here, the parameter p is related to the labels (p,p') of the minimal model as p = p/p'. 

The equation (|68|) has two linearly independent solutions having power-law behaviors 
at z — 0. These correspond to two conformal blocks in the S-channel, 



S*(z) 



1F1 



a + /3 + 7, a + (3' + 7 
+ a — a' 

a' + p + j, + + 7 
1 + a' — a 



(69) 



Another set of linear independent solutions having power law form at z — 1 corresponds 
to T-channel conformal blocks: 



T 1 (z)(z) = z a (l-zy 2 F X ( 



T 2 (z)(z) 



z a (l - z) 



2^1 



a + /3 + 7, a + /3' + 7 
1 + 7-7' 

a + P + j 1 , a + P' + j 1 
1 + 7' - 7 



1 - z 



(70) 



Since second-order differential equation 
S2 are expressed in terms of T\ and T 2 as 



has two linear independent solutions, S± and 



S 1 {z) = An(z) + BT 2 {z) , 
S 2 (z) = CT x {z) + DT 2 (z) . 



(71) 



Using the relations between solutions of Riemann equations (or just using the formulae 
for the analytic continuation for the hypergeometric functions), we obtain the expressions 
for the ^-independent coefficients 



A- 
B 
C 
D 



r(l + a-a / ) r (Y-7) 
T(a + p' + i)T(a + p + i) ' 

r(l + a - a') r (7 - 7') 
T{a + P + 1 )T{a + p' + 1 ) ' 

r(l + ct - a)Y{i - 7) 
T(a' + p' + 7 ')r(a' + /? + 7') ' 

r(i + a'-a) r( 7 -y) 

r(a' + /3 + 7)r(a' + /5' + 7) ' 



In what follows, we need the ratio AB/CD. It can be expressed in terms of the j(x) 
function as 



AB 
CD 



7(1 + a - a') 7(a' + P + Y)7(«' + /? + 7h(«' + + 7) 



7(1 + a' — a) 



j(a + P + 7) 
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For the conformal blocks in model Ai p , p >, the following expression holds for the operators 

AB _ i{n lP ) 7 (2 - (ni + n 2 + w 3 )f ) 7((~^i + n 2 + w 3 )f) 
CD 7(2 -n x p) 7((n 1 -n 2 + n 3 )^) 7((ni + n 2 - n 3 )|) 

The crossing symmetry condition [2] (see fig. H} leads to the equation for the correlation 

"1 "3 M n 3 

\ v\+\ / \ n, -1 / 

\ S2 




Figure 4: Conformal blocks. 

functions 

C(2,m,ni + l)C(ni + 1, "-2, ^3) 1 5*1 1 2 + C(2,ni,m - l)C(m - l,n 2) n3)|5 , 2| 2 = 
= C(2,n 2 ,n 2 + l)C(n 2 + 1, n 1; ng)!^] 2 + C(2,n 2 ,n 2 - l)C(n 2 - 1, n x , n 3 )|T 2 | 2 , 

where C(n 1; n 2 , n 3 ) denotes the structure constants C™- 5 of the operator algebra of fields 
In particular, this gives 

C(2, ni, ni + l)C(m + 1, "2, ras)AB + C(2, ni, m - l)C(ni - 1, n 2 , %)CD = . 
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We consider the case n 3 = ri\ = n and n 2 = 2. If A , B , . . . are the corresponding 
elements of the monodromy matrix, then we find 



C 2 (2, n,n + 1)A B + C 2 (2, n, n - 1)C D = . 



Therefore, 



C(2,n,n+ 1) 
C(2,n,n- 1) 



We can now write the general difference equation for the structure constants in the form 

C(ni + l,n 2 ,n 3 ) 



C(ni - l,n 2 ,n 3 ) 



A B i l 2 CD 

x 



C A, 



Rewriting everything via the variables rij, we obtain the functional equation for structure 
constants in the $i in subalgebra of M. v ^> 

C( ni + 2,n 2} n 3 ) _ n(2-(n 1 + l)p) 1 (2-(n 1 + 2)p) 



x 



C(ni,n 2 ,n 3 ) I jfap) j^m + l)p) 

7 ((m - n 2 + n 3 + l)f) 7 ((n a + n 2 - n 3 + 1)§) 



7 ((2 - m + n 2 + n 3 + l)f) 7((-ni + n 2 + n 3 - 1)0 ' 
We seek the solution in the ansatz 

J (hp) 



(72) 



C(ni, n 2 , n 3 ) = N(rii) J [ 



L 7((m - k)p)j((n 2 - k)p)^{2 - (n 3 + fc)p) 



where 



n = (ni + n 2 - n 3 - l)/2 . 
Substituting this expression in (J72J), we obtain the constraint for iV(ni) 

which can be easily resolved as 

N{n)= \ \{ 1 {kp) 1 {2-{k + l)p) 



>n-l 



(73) 



yk=l 
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Taking the symmetry A, < > Aj into account, we find the expression for the structure 

constants 

_ N( ni )N(n 2 ) " 7M , . 

lni ' n2 ' n3j N(n 3 ) ll 7 ((n 1 -%)7((n 2 -%) 7 (2-(n 3 + %) ' 1 J 

We briefly comment on the admissible triples (n 1; n 2 , n 3 ). The procedure used in the 
derivation assumes that the field 3>i, n3 must be related to $i, ni and $i, n2 by fusion 
rules [2J. Otherwise, the proposed structure constants should be reconstructed by sym- 
metry. 



B Integrals for the first corrections 

The integrals from the correlation functions that appear at the first order correction 
can be represented in the form of double integrals over a plane because of the equation 

0123 Em 



71" - 



j(a + 1)t(c+ 1) 



d 2 x\x\ 2 ^ a+d+c+1) \l-x 



2c 



+7T 



7(a + c + 2) 
7(6+1)7(0 + 0+1) , , 2i/ 



d xlxr II - x| 



7(a + 6 + c + 2) 
d 2 x / d 2 2/|^i 2a |l - xl 2b |y| 2d |l - 2/| 2e |a: - 2/ 



2F1 
2F1 



x 



a + 1, -6 
a + c + 2 

—a — 6 — c — 1, — c x 2 
—a — c 



x 



For practical computations it is convenient to decompose this integral into a sum 
of holomorphic and antiholomorphic parts using the way proposed in works (THl |Hj (see 
ITS] for alternative methods). Starting with the integral 



J 




d 2 xd 2 y \x\ 2a \l - x\ M \y\ M \l - y\ le \x - y 



1 2a 



26 1 .|2d| 



1 2e I 



1 2c 



(75) 



performing Wick rotation 

x 2 = ix (l - 2ie) , y 2 = iy {l - 2ie) 
and introducing the new variables 



X = X\ + Xq 

y = yi + yo, 



X — X\ — Xq j 

v = vi - Vo , 
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we can easily obtain 



J = - \jj dxdy x a (x - \fy\y - l) e (x - y) c 

x dxdy(x + ie(x-x)) a (x-l + ie(x-x)) b (y + ie(y-y)) d ( 76 ) 



x (y - 1 + ie(y - y)) e (x -y + ie((x - x) - (y - y))f . 
Integral (|75jl has the branch points 

x = —iex , 

x = l-ie{x - 1) , _ (77) 

_ _ . , v 1/ = 1 -ie(y - 1) • 
x = y- ie{x - y) , 

By deforming integration contours is easy to demonstrate that non-trivial contributions 
to (|7Sjl come only from two domains of integration, namely, 

{0 < x < 1 , < y < a;} and {0 < x < 1 , x < y < 1} . (78) 

One observes that the integral J is represented as sum 

J = I (a, b, d, e, c) + I(d, e, a, 6, c) , (79) 

where 

1 X 

if, r . ^ 



/(a, 6, d, e, c) = -- / dx / x a (x - 1) V (y - l) e (x - y) c x 







xy dx J dyx a (x-l)Y(y-l) e (x-y) c , 
Ci c 2 



1 1 

1 / , / , „ -\b„,d/ 



J(d, e,a } b } c) = -- I dx I dyx a (x - l)Y{v ~ l) e (^ - Z/) c x 



x 

\b-d/- 



x j dy j dxx a (x-l) b y d (y-l) e (x-y) c 

C4 C3 
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Figure 5: The contours C\, C 2 , C 3 , C 4 . 



We introduce notations I (a, b, d, e, c), I(d, e, a, b, c) for the integrals above since that two 
integrals are related by transform a — > d, b — > e as it can be easily checked by change of 
variables. The contours of integrations are shown in Fig. Let us now transform the 
contour integrals I(a, b, d, e, c) into standard type integrals. This is easily provided by 
applying the technique of analytical continuations and manipulation with the contours 
described in . We obtain the following expression 

I(a, b, d, e, c) = sin nbK(a, b, d, e, c) [sin TreL(a, b, d, e, c) + sin 7r(e + c)L(d, e, a, b, c)) , (80) 

where 

i i 

K = J dx J dyx a+d+c+1 {\ -z)V(l -y) c (l ~xy) e , 



1 1 

L = J dy J dx x~ a - b - c - 2 (l - x) c y- a - b - d - e - c - 3 < 



Now one uses standard formulae for the integral representations of the hypegeometric 



51) 



\i- y y(i- xy y 
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functions 
i 



dy y p {l - yf{l - xy) r = B(p + 1, q + 1) 2 F t 



P + l 



p + q + 2 



x 



J x s {l-xf 2 FiQ' 6 |x)rfx = fl(s + l,t + l) 3^2 ( 



6, s + 1 



s + t + 2, c 



to find that the integrals K and L admit another representation in terms of the higher 
hypergeometric functions. Here B(x,y) = T(x)T(y)/T(x + y) is a beta function. Namely, 



K = B(d + 1, c + l)5(a + rf + c + 2,6 + l) 

— e, a + <i + c + 2, d+1 

a + d + c + 3 + 6, d + 2 + c 



X 3^2 



(82) 



L = B(-a - b-d-e-c-2,e + l)B(-a -b- c- l,c + 1) 

—6, —a — 6 — d — e — c — 2, —a — b — c — 1 

—a — b — d — c — 1, —a — b 



X 3^2 



1 . 



Equations (JTUJ) , (IHOJ) and (JS2|) give a form convenient for numerical studies. We note 
that the functions 3F2 at unity satisfy specific identities (37] • This simplifies the answer 
and leads to 
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